Storing entanglement of nuclear spins via Uhrig Dynamical Decoupling 
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Stroboscopic spin flips have already been shown to prolong the coherence times of quantum sys- 
tems under noisy environments. Uhrig's dynamical decoupling scheme provides an optimal sequence 
for a quantum system interacting with a dephasing bath. Several experimental demonstrations have 
already verified the efficiency of such dynamical decoupling schemes in preserving single qubit co- 
herences. In this work we describe the experimental study of Uhrig's dynamical decoupling in 
preserving two-qubit entangled states using an ensemble of spin- 1/2 nuclear pairs in solution state. 
We find that the performance of odd-order Uhrig sequences in preserving entanglement is superior 
to both even-order Uhrig sequences and periodic spin-flip sequences. We also find that there exists 
an optimal order of the Uhrig sequence using which a singlet state can be stored at high correlation 
for about 30 seconds. 
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I. INTRODUCTION 

Harnessing the quantum properties of physical sys- 
tems has several potential applications, particularly in in- 
formation processing, secure data communications, and 
quantum simulators It is believed that such quan- 
tum devices may play an important role in future tech- 
nology [2]. But their physical realization is challenging 
mainly because of decoherence - the decay of the coherent 
states due to interaction with the surrounding environ- 
ment 01 ■ Therefore it is important to minimize the 
effects of decoherence using suitable perturbation on the 
quantum system [sj. One technique known as 'dynam- 
ical decoupling' involves protecting the quantum state 
from decoherence by driving the system in a systematic 
manner such that the effective interactions with the envi- 
ronment at different instants of time cancel one another. 

Preserving nuclear spin coherences by spin flips at reg- 
ular intervals was long been known in nuclear magnetic 
resonance (NMR) as the famous Carr-Purcell-Meiboom- 
Gill (CPMG) sequence [1, 0. The CPMG sequence is 
widely used in NMR to measure the transverse relaxation 
time constants in the presence of spatial inhomogeneity of 
the static magnetic field and temporal fluctuations in the 
local fields arising due to molecular motion [8,]. The se- 
quence involves a set of N tt pulses uniformly distributed 
in a duration [0, T] at time instants {^1,^2, • • • , ^at}- As- 



suming instantaneous tt pulses, 
in j 
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time instant is linear 



(1) 



Of course, in practice the tt pulses do have finite dura- 
tion owing to the limited power of electromagnetic ir- 



radiation generated by a given hardware. Further, the 
constant time period between these spin flips should ide- 
ally be shorter than the correlation time of the spin-bath 
interaction. Even this delay is limited by the maximum 
duty-cycle that is allowed for the hardware. Dynamical 
decoupling with such bounded controls have also been 
suggested [9|-[T2|. For instance Hao et al. have been able 
to calculate, using a particular type of atomic systems, 
the maximum delay between spin-flips in order to effi- 
ciently suppress decoherence due to a bath with a finite 
cut-off frequency [l^ . By studying the efficiency of the 
decoupling as a function of the CPMG period often it 
is possible to extract valuable informations about molec- 
ular dynamics and such studies are broadly categorized 
under 'CPMG dispersion' experiments 

Uhrig generalized the CPMG sequence by considering 
an optimal distribution {ti, t2, ■ ■ ■ ,1^} of N spin flips in 
a given duration [0, T] of time that provides most effi- 
cient dynamical decoupling J^]. Using a simple dephas- 
ing model, Uhrig proved that the time instants should 
vary as a squared sine bell: 
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UDD works well in systems having a high-frequency 
dominated bath with a sharp cutoff |16l - [l8j . On the 
other hand when the spectral density of the bath has a 
soft cutoff (such as a broad Gaussian or Lorentzian) , the 
CPMG seq uence was found to outperform the UDD se- 
quence |19l - [25| . Suter and co-workers have studied these 
different regimes and arrived at optimal conditions for 
the dynamical decoupling (26j . 

Recently Agarwal has shown using theoretical and nu- 
merical calculations that even entangled states of two- 
spin systems can be stored more efficiently using UDD 
[27| . Since entangled states play a central role in QIP, 
teleportation, data encryption, and so on, saving entan- 
'^Sshfflgfe.iis^aiiQafeEai.tm ejjilcysical realization of quantum de- 
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FIG. 1: (Color online) The NMR spectrum and the molec- 
ular structure of 5-chlorothiophene-2-carbonitrile. 



vices More recently dynamical decoupling on an 

electron-nuclear spin-pair in a solid state system has been 
shown to prolong the pseudoentanglement lifetime by two 
orders of magnitude [28| . 

While much of the experimental effort has been on test- 
ing the loss of coherence due to T2 processes, ours is one 
of the first experiments where we study not only the loss 
of coherences, but also the loss of entanglement due to 
both Ti and T2 processes. Though newer sequences have 
been suggested to decouple both of these processes, these 
are yet to be studied experimentally 29,, ^] . Since we 
have developed experimental techniques where we can 
prepare Bell states with high fidelity and characterize 
these states with high precision [sil, [s^l , we explore the 
utility of different dynamical decoupling sequences on 
systems wherein both Ti and T2 relaxations are signif- 
icant. In the following we describe preparation of Bell 
states using a two-qubit NMR system and then describe 
the comparative performances of CPMG and UDD on 
such entangled states. 



II. PREPARATION AND STORAGE OF 
ENTANGLEMENT 

We study storage of entanglement by dynamic de- 
coupling on a pair of spin- 1/2 nuclei using liquid state 
NMR techniques. The sample consisted of 5 mg of 
5-chlorothiophene-2-carbonitrile dissolved in 0.75 ml of 
dimethyl sulphoxide (see Figure [1]). The two protons of 
the solute molecule differ in the Larmor frequency by 
Ai/ = 270.4 Hz and have an indirect spin-spin coupling 
constant of J = 4.1 Hz. The T2 relaxation time constants 
for the two protons are about 2.3 s and the Ti relaxation 
time constants are about 6.3 s. All the experiments are 
carried out in a Bruker 500 MHz NMR spectrometer at 
an ambient temperature of 300 K. 



Pseudopure States 

An ensemble of nuclear spin-systems in thermal equi- 
librium at room temperature T exists in a mixed state at 
high spin temperature, since the Zeeman energy gaps /S.E 
are much smaller than the Boltzman energy kT. There- 
fore the general state of an ensemble of spin- 1/2 nuclear 
pairs can be described by the density matrix 



-1 



(3) 



where the identity part is usually ignored and the trace- 
less matrix is called the deviation density matrix. The 
dimensionless quantity e ^ IS.E /kT and have magnitudes 



~ 10~^ for protons in currently available magnetic fields 
at ordinary room temperatures. Using a combination of 
unitary and non-unitary processes, it is possible to trans- 
form the above state to a pseudopure j^A) state given by 



Ppps = |(l-e)l + e|V')(^l 



(4) 



[33| . Here e is a measure of the magnetization retained 
in the pseudopure state. All the kcts in the rest of the 
article must be understood as pseudopure states. 



Preparation of pseudoentangled States 

High fidelity entangled states are prepared via long 
lived singlet states in a procedure described by Soumya 
et al [31, 32]. The method is briefly explained in the 
following for completeness. 

The Hamiltonian for an ensemble of spin- 1/2 nuclear 
pairs of same isotope, in the RF interaction frame, can 
be expressed as 



cff 



Aj/ 



/2 + JI^ ■ /2 



1^12 



Jl,2 



(5) 



Here h/S.v corresponds to the difference in the Zeeman 
energy gaps of the two spins (due to the chemical shift 
difference), corresponds to the spin angular momen- 
tum of fc"^ spin, ,/ is the indirect spin-spin coupling, and 
1^12 is the amplitude of circularly polarized component of 
the radio frequency field, which is at resonance with the 
mean frequency of the two spins. 

In the limiting case of lS.v — ^ 0, the system is said to 
have magnetic equivalence, and the singlet state |5o) = 
(|Gl)-|10))/\/2, and the triplet states \Ti) = |00), \Tq) = 
(|01) -I- |10))/%/2, and |T_i) = |11) form an orthonormal 
eigenbasis of the resulting internal Hamiltonian Q 



(6) 



The equivalence Hamiltonian may be realized by sup- 
pressing the chemical shift Ai/ using a radio frequency 
'spin-lock' m. 

The long-lived nature of singlet states under the equiv- 
alence Hamiltonian can be used to prepare high-fidelity 
Bell states. The experiment involves preparing an inco- 
herent mixture of singlet and triplet states 



P(0) 



|^o>(^o|-|To)(To| 



(7) 



from the equilibrium state by using the pulse se- 
quence shown in Figure [2] [3J] . After preparing p(0), the 
Hamiltonian 7^°^ is imposed using RF spin-lock. During 
the spin-lock \Tq) state rapidly equilibrates with the other 
triplet states. On the other hand, the decay constant 
of singlet state j^o) during the spin- lock is much longer 
than the spin-lattice relaxation time constant (Ti) (and 
hence the singlet state is known as a long-lived state) 
[ssl . [36| . The long-life is attributed to the fact that the 



3 



n/2. 



TT/2. 



1 4J I 4J 2Av 


|_L 




|So>*|i|J,> 
|So>*l<P.> 


/ s 




4Av 


Spin-lock 


UDD-N 


Tomography 


» i (optional) 







i; + 1/ 



|S„>,|T„> |S„> 



FIG. 2: (Color online) NMR pulse sequence to study dynami- 
cal decoupling on Bell states. An incoherent mixture of singlet 
and triplet states is prepared which under spin-lock purifies 
to singlet state. The resulting singlet state can be converted 
to other Bell states. Then dynamical decoupling sequence can 
be applied and the performance of the sequence can be stud- 
ied by characterizing the residual state using density matrix 
tomography. 



anti-symmetric singlet state can not be converted to sym- 
metric triplet states via symmetry conserving relaxation 
mechanisms such as the intra-molecular dipolar interac- 
tions d^-llii- 

As a result of the long life-times under , the atten- 
uated correlation of the singlet state 



(So) 



trace[|^o)(^o| -pft)] 
i/trace [p'^(t)] 



(8) 



with p{t) improves over time t, reaches a maximum value 
and ultimately decays due to the symmetry breaking in- 
teractions [3l[. The spin- lock imposing the Hamiltonian 
"H^q can be turned off when the maximum correlation 
(So) is reached, and the high fidelity singlet state can be 
used for further studies. 

Other Bell states can be obtained easily from the sin- 
glet state: 



^ |V'+) = (|01) + |10))/V2, 
^ !</)_) = (100) -|11))/V2, 
^ |0+) = (|OO) + |11))/V2. (9) 



l^o) 
l^o) 
l^o) 



The z-rotation in the above propagators can be imple- 
mented by using chemical shift evolution for a period of 
l/{2Av), and qubit selective x- rotation can be imple- 
mented by using radio frequency pulses [32[ . 

Details of dynamical decoupling on the Bell states will 
be described in the next sub-section. In order to inves- 
tigate the decoupling performance one needs to quantify 
the decay of Bell states with decoupling duration. The 
Bell states by themselves are inaccessible to macroscopic 
observables, but can indirectly be detected transforming 
to observable single quantum coherences [SSj, i36j . Alter- 
natively, a more detailed and quantitative analysis of Bell 
states may be carried out using density matrix tomogra- 
phy [sij . After measuring the density matrix, attenuated 
correlation can be evaluated for any Bell state using ex- 
pressions similar to (|8]). 

In the following we describe experimental implemen- 
tations of dynamical decoupling after preparing such en- 
tangled states. 
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FIG. 3: (Color online) Pulse sequences for various orders of 
Uhrig Dynamical Decoupling. Note that both UDD-1 and 
UDD-2 are equivalent to CPMG. The time instants are cal- 
culated according to the expression ([2]), with A'' being the 
order of UDD and the total period T = N x 4.0272 ms. 



Dynamical Decoupling 

As described earlier, the UDD scheme consists of a se- 
quence of spin flips placed at time instants given by the 
expression ([2]). Instead of applying the Uhrig's formula 
for the entire duration of decoupling, we apply the for- 
mula for a short time interval (T) consisting of a small 
number (N) of pulses and then repeat the sequence. Fig- 
ure [3] shows pulse sequences for various orders of Uhrig 
Dynamical Decoupling (we refer to an N-pulse UDD se- 
quence as UDD-N). Note that UDD-1 (and UDD-2) are 
equivalent to CPMG sequence, in which repeating seg- 
ment consists of [rcpMG — tt — tcpmg]- In our experi- 
ments, TcpMG was set to 2 ms and the duration of the 
TT pulse was 27.2 fis. The total duration of UDD-N was 
set to T = 7V(2rcpMG + Ttt), such that for an extended 
period of time, the total number of tt pulses remains same 
irrespective of the order of UDD. Only the distribution 
of TT pulses varies according to the order of UDD. For 
example in one second of decoupling, there will be about 
250 TT pulses in all UDD-N. Our investigation thus helps 
in studying the efficiency of decoupling over a fixed du- 
ration of time for a given number of n pulses dispersed 
according to different orders of UDD. 

Now we describe the performances of UDD-N on the 
singlet state which was prepared as explained before (see 
Figure [2]). After applying UDD-N for a fixed duration 
of time, we carried out density matrix tomography and 
evaluated the correlation of the preserved state with the- 
oretical singlet density matrix. The correlations for var- 
ious orders of UDD are displayed in Figure HI As can be 
seen from the figure, the singlet state can be preserved for 
longer durations by UDD-1 (CPMG) than no-decoupling. 
It is also clear that all even-order UDD sequences result 
in significant fiuctuations in the correlation of the singlet 
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FIG. 4: (Color online) Experimental correlations (circles) of 
singlet state as a function of decoupling duration of various 
orders of UDD. Also shown in the top-left figure is the corre- 
lation decay under no dynamical decoupling (squares). 



state. However, the odd order UDD preserve the singlet 
state for tens of seconds. For example, the correlation of 
the singlet state under UDD- 7 at all the sampled time 
points till 20 seconds is above 0.96. This rather surpris- 
ing even-odd behavior is likely due to the differences in 
the performances of the even and odd ordered sequences 
against the spatial inhmongeneity of the RF pulses. 

One way to quantify the efficiency of dynamical de- 
coupling under various orders of UDD in figure lU is by 
counting the number of time instants in which the cor- 
relation of the preserved state exceeds a given threshold. 
The bar plot in Figure [S] compares the number of time 
instants during decoupling under various orders of UDD 
in which the correlation of the singlet state exceeded 0.9. 
ft can be seen that there exists an optimal order of UDD 
(for a given tcpmg and r^), which performs the most ef- 
ficient decoupling. The optimality may be because of the 
finite width of the tt pulse. In a CPMG sequence the tt 
pulses are uniformly dispersed, while in Uhrig sequence 
the TT pulses are more crowded at the terminals (begin- 
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FIG. 5: (Color online) The number of time instants at which 
the correlation exceeded 0.9 for various orders of UDD-N. 
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FIG. 6: (Color online) The decay of the singlet spin-order 
measured by converting it into observable single quantum 
magnetizations. The decay was studied under CPMG se- 
quence (squares) as well as under Uhrig sequence (filled cir- 
cles). The dashed and the solid line correspond to the expo- 
nential fits for CPMG and UDD-7 data points respectively. 



ning and ending) of the sequence. For example, if there 
are too many tt pulses, Uhrig's formula will lead to an 
overlap of pulses. Experimentally, the overcrowding of 
TT pulses may also lead to RF heating of the sample and 
the probe. Thus the performance of the UDD sequence 
does not grow indefinitely with the order of the sequence, 
but instead will fall beyond a certain order. In our ex- 
perimental setting, we find that UDD-7 is the optimal 
sequence for storing the singlet state. There are recent 
suggestions for decoupling using finite pulses, however 
these are yet to be studied experimentally [i^, SH . 



It can be noticed that the attenuated correlation (ex- 
pression ([8])) is insensitive to the decay of the overall 
magnetization (e in (|3])), but simply measures the over- 
lap between pA and the theoretical density matrix 
An alternate method is to monitor the decay of magne- 
tization (i.e., e) under dynamical decoupling. As already 
mentioned, singlet state itself can not be measured di- 
rectly, but can be converted to observable magnetiza- 
tion by using a chemical shift evolution for a duration 
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FIG. 7: (Color online) Experimental correlations of the prod- 
uct state and various Bell states as a function of duration 
under (i) no decoupling (open squares), (ii) CPMG sequence 
(filled circles), and (iii) UDD-7 (open circles). 

l/(4Ai/) followed by a (f pulse. Intensity of the re- 
sulting signal as a function of tlie duration of dynamical 
decoupling is shown in Figure [5] As can be seen, UDD-7 
is no better than CPMG in preserving the overall spin- 
order. In fact the decay constant for CPMG and UDD-7 
are 6.1 s and 5.9 s respectively. 



Now we compare the efficiency of the optimal sequence 
UDD-7 with UDD-1 (CPMG) for preserving product 
state (cr^ + (7^) and other Bell states. Figure [7] shows the 
variation of correlation of product states and the Bell 



states as a function of the decoupling duration. Here, 
after preparing each of the initial state, the dynamical 
decoupling was applied for a fixed duration of time. To 
monitor the correlation, we have carried out the density 
matrix tomography as described earlier [3l|. In the case 
of no decoupling, we observe a rapid decay of the corre- 
lation. The UDD-1 (CPMG) sequence shows some im- 
provement in the storage time. However, UDD-7 clearly 
exhibits much longer storage times than the CPMG se- 
quence. The superior performance of UDD-7 on the sin- 
glet state compared to other Bell states is presumably 
because of its antisymmetric property described in sec- 
tion II. 



III. CONCLUSIONS 

We have studied the efficiencies of CPMG as well as 
Uhrig dynamical decoupling sequences on 2-qubit Bell 
states both in terms of magnetization as well as in terms 
of correlation decay. While the Uhrig sequence is no bet- 
ter than CPMG sequence in terms of preserving the over- 
all magnetization (or spin order), it clearly outperforms 
the CPMG sequence in preserving the correlation of the 
entangled as well as non-entangled states. We summa- 
rize three important features: (i) the even-order UDD 
sequences result is fluctuations in correlations, (ii) the 
odd-order UDD sequences out-perform the CPMG se- 
quence, and (iii) there exists an optimal length for the 
odd order UDD sequence which exhibits most efficient 
decoupling. In our case, UDD-7 of 28.2 ms duration ap- 
peared to outperform all other sequences of both lower 
and higher orders. We are carrying out investigations 
into the effects of other experimental issues like RF in- 
homogeneity, resonance off-set, errors in calibration of 
pulse angle etc. These considerations may help in the 
theoretical and practical understanding of the optimal 
decoupling schemes. 
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